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Abstract 

The cross sections for charm and bottom quark production in the 
threshold region are discussed. We consider the effects of an all order 
resummation of initial state soft-plus-virtual gluon radiation on the 
total cross sections compared to the order a% results. 
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Heavy quark production has long been a topic of interest, both experi- 
mentally and theoretically. Early measurements of the total cc production 
cross section at \/S < 63 GeV suggested that the calculated Born (LO) cross 
section underpredicted the data by a factor of two to three [1], |J, known as 
the K factor after a similar situation in Drell-Yan production. In general, 

K _ cr d ata(AB ->■ QQ) 

6XP f7 t heory(^4-B ~> QQ) 

where Q is the heavy quark (c, b, or t) and cr t heory is calculated to fixed 
order in the running coupling constant a s . The next-to-leading order (NLO) 
corrections to the Born cross section have been calculated [|, |], |5| and an 
analogous theoretical K factor is defined by the ratio of the NLO to the LO 
cross sections, 

a^°(AB -> QQ) 
th ~ a(°)(AB -> QQ) ' K } 

where is the LO cross section and cr NLO is the sum of the LO and the 
exact 0(a s ) correction, cx NLO = + | ex act- 

The heavy quark production cross section is calculated in QCD by as- 
suming the validity of the factorization theorem || and expanding the con- 
tributions to the amplitude in powers of the coupling constant a s (/i 2 ). The 
hadronic production cross section at hadronic center of mass energy y/S and 
to order a{ k ^ is 



ij J — Jt x x 



where f^x,^ 2 ) are the scale-dependent parton densities of hadron h eval- 
uated at the scale /i 2 and a^j is the k th order partonic cross section for 
ij — > QQX. The total cross section up to order k is J2k c^ k '(S, m 2 ). The 
numerical results for "lighter" heavy quark production depend on the choice 
of the parton densities, involving the mass factorization scale fip, the run- 
ning coupling constant, evaluated at the renormalization scale hr, and the 
heavy quark mass m. Usually the renormalization and factorization scales 
are assumed to be equal, //r = \if = \i. All these factors influence both K exp 
and K th . 

At LO the cross section is very sensitive to the mass and scale parameters. 
However even including the NLO corrections cannot completely fix the cross 
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section. The sensitivity to even higher terms in the QCD expansion is often 
demonstrated by varying /x between m/2 and 2m. This may not be very 
meaningful, especially for charm, since a variation of an order of magnitude 
or more is observed. It is therefore not clear that the next-to-next-to-leading 
order (NNLO) corrections are not at least as large as the NLO corrections, 
particularly when m <C y/S. Given these facts, it is impossible to make more 
precise predictions in the absence of a NNLO calculation. 

Recently two groups have attempted to use the NLO calculations to make 
more definitive statments about charm production J7|, || |§. The NLO calcu- 
lations were compared to charm production data in pp and ir~p interactions 
10| , PI - In the first approach [^, ||, the data were used to fix m c and \i 



by requiring ~ 1 in an attempt to place bounds on charm production 

at nuclear colliders. Two r ecent parton densities^, MRS D-' |3|, H and 
GRV HO [|T5|, were used in the comparison. Reasonable agreement was 
found for m c = 1.2 GeV/c 2 , // = 2m c with MRS D-' [|IJ, [Lj and for m c = 1.3 



GeV/c , fj, = m c with GRV HO [Jy| [16[ although both results tend to un- 
derestimate the total cc production cross section, crjg*, with K^P ~ 1.1 — 2. 
In the range of the parameter space defined by m c , fiR and /if, can 
be reduced to unity. However, it is questionable if the mass and scale values 
needed for K^!± ~ 1 are consistent with a perturbative treatment and with 
the defined limits of the parton density distributions. In another approach, 
calculations using a charm mass of m c = 1.5 GeV/c 2 produced results com- 
patible with the data although with some essential caveats: /ip and hr were 
varied independently and out-of-date parton distributions fit with several val- 
ues of Aqcd were used 0. Decreasing with respect to fiF and increasing 
Aqcd result in significantly larger cross sections for a given m c . Addition- 
ally, different parton densities were used in the calculations of a l °^ and high 
energy b production. Neither approach is fully satisfactory: either an un- 
comfortably small charm quark mass is needed or the parameters used to 
describe low energy production are incompatible with those used at collider 
energies. 

Although a complete calculation of still higher order terms is not possible 
for all values of S and m, improvements may be made in specific kinematical 
regions. Investigations have shown that near threshold there can be large 



2 A11 available parton densities for the nucleon and the pion can be found in PDFLIB 
HI. 
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logarithms in the perturbative expansion which must be resummed to make 
more reliable theoretical predictions. These large logarithms arise from an 
imperfect cancellation of the soft-plus-virtual (S+V) terms. In [17 an ap- 
proximation was given for the S+V gluon contributions and the analogy with 
the Drell-Yan process, studied in |L8|, was exploited to resum these to 
all orders of perturbation theory. The same resummation procedure was also 
applied to cr*-* and inclusive top quark distributions in pp collisions at the 
Fermilab Tevatron 
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For top production in pp collisions, qq an- 
nihilation is the dominant process, fortunate since the exponentiation of the 
S+V terms |T7J is better understood because the simple color structure has 
a close correspondence to the Drell-Yan studies |L8L figfl . 



The resummation of the leading S+V terms JTTJ modifies eq. (3) so that 



a 



5 (S,m 2 



ij Jtq Jt X X 



(tS, m 2 



where 



a 



ij 



! (rS, m 2 



-2ms 1 / 2 



ds 4 f( 



s 4 m 



s 4 , m 2 ) 



d,S4, 



(4) 



(5) 



The integration variable s 4 , an invariant which measures the four-momentum 
carried away by the final-state gluon in the process i(ki) + j(k 2 ) — ► Q(pi) + 
Q{pz) + g(k 3 ), is defined as s 4 = s + t + u — 2m 2 where s = (k± + k 2 ) 2 , 
t = (ki — pi) 2 , and u = (hi — P2) 2 ■ The cross section 0$ (s, s 4 , m 2 ) is 

the angle-averaged Born cross section and a\f (s,0,m 2 ) = a^(s,m 2 ). The 
function /(s 4 /m 2 , m 2 / n 2 ) is 



( S4 


m 2 \ 
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= exp 




\m 2 






7T 


\m z J 


771 . 



[s 4 m 2 } 71 

— — -exp {-r) lE , 6 
r(l + r/) 

where 7^ is the Euler constant, 77 = (8CV,-/ (3q) ln(l+(/? «s(/^ 2 )/47r) ln(m 2 / /i 2 )), 
Po = 11 — 2n//3 for SU(3), and A and a s are scheme dependent. In e.g. the 
MS scheme, A = 2 and a s ((s 4 /m 2 ), m 2 ) = a s (s 2 / 3 m 2 ^). The integral over 
s 4 has been split into two parts and the integral over the region < s 4 < s , 
where the running coupling becomes infinite as s 4 — > 0, is assumed to be 
negligible. This condition is satisfied if So/m 2 <C 1. The lower limit on the r 
integral is 



[m+ (m 2 + s ) 1/2 ] 2 
S 



(7) 
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where sq = m 2 ^^/ /i 2 ) 3 / 2 in the MS scheme and sq = m 2 {jj^/ ' /j 2 ) in the 



DIS scheme. The cutoff parameter /ig, introduced in [19, monitors the 



sensitivity of the cross section to nonperturbative higher-twist effects. If cr|J s 
is strongly dependent on /x , a precise determination of the cross section 
requires full knowledge of the nonperturbative contributions. As may be 
expected, the resummed corrections diverge for small /x [0- 

The resummation technique has recently been applied to bb production 
in pp interactions @ at HERA-B (^/S = 39.2 GeV) ||, [8J where gluon 
fusion is the dominant channel. The gg channel is not as amenable to the 
resummation technique because each of the Born diagrams has a different 
color structure^]. Therefore it is important to test the model in a regime 
where data already exists. Charm and bottom production data from pp and 
7i~p interactions are available at energies where the resummation technique 
can be applied. An advantage of the pion beam is the relative importance of 
the qq channel in cc and bb production. However charm production must be 
treated with some care. The data is available for v^S* > 15 GeV where the 
expansion parameter, a s (m 2 ) \n(y/S/m), is not small. However, if the model 
is reliable for charm production, a better understanding of cc production at 
the lower fixed-target energies may be reached. 

Since we wish to compare our results with data taken using both pion 
and proton beams, we are rather constrained in our choice of parton dis- 
tribution functions. The only consistent NLO evaluation of the pion and 



proton parton densities has been done by GRV |L5|, |1(|. Because the GRV 
HO parameterization is in the MS scheme, we are forced to use this scheme 
for both the qq and gg channels. Thus the value of /io needed for stability 
of the resummed result in the qq channel is likely to be larger than found 
previously ]17|, |(], |2l], p3|[ . Our results for the exact, approximate, and re- 



summed hadronic cross sections will be calculated using these distributions 
along with the two-loop uncorrected running coupling constant^. This set 
has the additional advantage of a rather small initial scale so that we can 
use \i = m c in our calculations. We therefore take m c = fi = 1.5 GeV/c 2 and 
mj = /i = 4.75 GeV/c 2 , along with the GRV HO parton densities and the MS 
scheme for our principle results. We will also show the range of resummed 

3 Quark-gluon scattering, producing a final state quark or antiquark, cannot be re- 
summed by this method since there is no equivalent Born term. 

4 The difference between the uncorrected running coupling constant and the corrected 
value given by PDFLIB O is small, « 3-4%. 
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cross sections by varying all these parameters, including the parton densities. 
In Table 1 we give the value of A n/ , \i and a s for each value of m and set of 
parton densities we consider. The choice of these particular parameters will 
be explained later. 

We first review the NLO contributions to cc and bb production || £|, |]. 
In Fig. 1 we show the relative contributions of the gg, qq, and the (q + q)g 
channels as functions of v^S* for cc and bb production in n~p and pp interac- 
tions. It is clear that cc production is almost entirely dominated by gluon 
fusion. In ir~p interactions at \/S ~ 10 GeV the gg channel accounts for half 
the cross section. However, the gluon contribution increases to more than 
80% soon afterward and the gg channel in pp production is always ~ 90% of 
the total cross section. Note that the quark-gluon scattering contribution is 
negative at low energies but changes sign at \/S ~ 40 GeV. The qq annhila- 
tion contribution is larger for bb production, particularly for 7r~p interactions 
where this channel is dominant until \/S ~ 40 GeV. This situation is similar 
to top production in pp collisions at the Fermilab Tevatron. The pp pro- 
duction cross section is again dominated by gluon fusion, between 70-80% 
of the total, somewhat less than the gg contribution to cc production at the 
same energies. Quark-gluon scattering always gives a negative contribution 
to the total bb cross section for the energies considered. Note that although 
we show the bb results for y/S = 10 GeV, the calculation is most reliable for 
VS > 20 GeV. 

The resulting theoretical K factors for the total cross sections reflect the 
uncertainties in the dominant channels. Thus for cc production, > 2, 
near the gluon K t h, over a large energy range. Likewise, K t ^ is smaller for 
n~p — > bb interactions because of the dominant qq contribution for v^S* < 
40 GeV. The qq channel is most amenable to resummation because even 
above threshold the NLO corrections are small. Above threshold in the 
gg channel, the NLO contribution increases as a function of \/S/2m and 
becomes constant at high energies while the Born contribution decreases to 
zero, resulting in large corrections. In bb production K t h is large for both 
channels near threshold. 

In Fig. 2 we examine the /xo dependence of the resummed cross section 
using our parameters for cc production at y/S = 15 GeV and bb production 
at y/S = 30 GeV. We also show, for comparison, the /Jo dependence of the 
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approximate NNLO cross section, 



a app^(0) + (J (l) |app+(T (2) |app; (g) 

where we have imposed the same S4 (34 > sq) phase space cut on <r app as on 
the resummed cross section. Here a^ l > | app and a^ 2 > | app denote the approx- 
imate first and second order corrections, respectively, where only soft gluon 
contributions are taken into account. Note that the phase space cut is also 
applied to a^°\ The effect of the resummation is apparent in the difference 
between the curves. At small fiQ, a ms diverges, signalling the onset of non- 



perturbative physics. In a study of bb production at HERA-B p3| , /i was 
chosen so that the resummed cross section remains somewhat larger than 
cr app in each channel, namely fi = 0.13m in the qq channel with the MRS 
D— ' DIS distributions and /io = 0.36m in the gg channel, quite similar to 
earlier choices for top quark production pOflPl- We choose the same ratios p3 



for the gg channel, so that /io ~ 0.35m. Reasonable convergence is seen for 
this value, even for cc production, although <r res remains larger than <r app for 
all values of Ho/m shown. Even though our calculation is in the MS scheme, 
convergence is found for /x ~ 0.15m in the qq channel, not significantly larger 
than that of the DIS scheme. 

In Fig. 3 we plot the resummed cross section with our chosen values of 
fiQ as a function of center of mass energy for cc and bb production at fixed- 
target energies. Since the exact 0{a^) results are known, we also show the 
perturbation theory improved cross sections defined by 

a imp = a res + a« | exact -a« | app , (9) 

to exploit the fact that | cxact is known and | app is included in cr rcs . 
The difference between <r res and a imp is larger in pp production, presumably 
because the qq approximation of <J^\ cxac t is better than the gg approxima- 
tion. We also show the NLO cross section calculated with the same mass and 
scale factors as <7 res and <j imp . The resummed and NNLO approximate cross 
sections were calculated with the cut s 4 > s while no cut was imposed on the 
NLO result. We can expect the resummation technique to work well in both 
the qq and gg channels when m/y/S < 0.1 because the NLO contribution 

5 An application of principal value resummation leads to very similar results for the top 
quark production cross section [Em . 
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is small compared to the LO term and the perturbation series should then 
converge. The effect of resummation is large for cc production, increasing 
the cross section by a factor of five or more relative to the NLO result. The 
agreement with the cc data from ix~p (I], ||, [Li] and pp [p], §, [UJ interactions 



is improved by the resummation. The agreement with the existing bb data 
from ir~p |27| interactions is also somewhat improved. In general, we see that 
Kl^ p is much smaller than for the same m, /i, and parton densities. 

However, for charm quark production, the technique is questionable above 
= 15 GeV. The expansion parameter is proportional to a s (m 2 ) \a{\fS/m) 
so that when 10 < v^S < 30 GeV, 0.53 < a s (m 2 c ) \n(^S/m c ) < 0.83 where 
a s (m 2 ) for the GRV HO distributions is found in Table 1. We have shown 
the cc results up to v^S* = 30 GeV even though the perturbative expansion 
no longer converges and resummation fails in the gg channel. This can be 
clearly seen in the faster increase of <r res and <r imp with energy compared to 
(x NLO for > 20 GeV in Fig. 3(a). 

We have also checked the resummation technique for a range of heavy 
quark masses, scales, and parton densities. The variation in the results is 
indicated by the dotted curves in Fig. 3. Outside the range indicated by 
the dotted curves, the resummation technique becomes questionable, par- 
ticularly for the bound given by the upper dotted curve. When examining 
bb production, we studied the mass and scale dependence over the range 
4.5 < mf, < 5 GeV/c 2 and m&/2 < < Im^. The bb results are more stable 
with respect to the variation, as indicated by the narrower bands in Fig. 
3(c) and (d). We varied the charm quark mass between 1.3 and 1.8 GeV/c 2 . 
Since the charm quark mass is relatively light, we study the scale dependence 
indirectly. Note also that the resummation technique cannot be applied to 
scales where [i < m c because a perturbative treatment is uncertain for mass 
scales of ~ 1 GeV/c 2 . 

We also studied the parton density dependence so that along with the 
GRV HO distributions we also used the GRV LO and the MRS D— ' proton 
distributions. Since the resummed cross section depends on it is in- 
structive to compare the results using LO parton densities with the one-loop 
coupling constant to the NLO parton densities with the two-loop running 



coupling. We also used the MRS D-' |J and SMRS P2 g3| pion densi- 
ties for MS 7c~p production. The low x behavior of the SMRS densities is 
different than the MRS D— ' densities. Presumably adjusting the behavior 
of the SMRS low x region to match that of the proton densities would shift 
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the large x contribution in the threshold region. Additionally, because Aqcd 
is different for the two sets we use the value obtained in the proton fit. Fi- 
nally for the MRS distributions the initial scale is larger than m c , Ql = 5 
(GeV/c) 2 . Therefore when we calculate cc production with these densities, 
we take fi = 2m c instead of a full scale variation while we use fi = m c for the 
GRV distributions. 

The larger value of a s at one-loop both increases the cross section with the 
GRV LO distributions as well as making the convergence of the perturbation 
series slower for both cc and bb production. In fact, since a s ~ 0.38 at 
one-loop for m c = 1.5 GeV/c 2 , the resummation technique appears to fail 
for the gg case even at this relatively low energy, y/S = 15 GeV. In general, 
larger values of fi makes the convergence of the perturbation series unreliable. 
When /i = 2m, s is reduced by a factor of four, resulting in a large a s in the 
exponent in eq. (6). Thus a much larger /i is needed for convergence, /i ~ 
0.5m c in the qq channel and /j,q ~ 0.7m c in the gg channel with the MRS D— ' 



distributions, much too large to be compatible with the earlier results [[17 



|20| , |23|p|. We also compared the MS and DIS schemes for the qq channel 
when appropriate. For both the GRV LO and MRS D— ' DIS distributions 
(applied to pp production only since the SMRS pion distributions are only 
available in the MS scheme) we found convergence for a smaller value of /z 
and a correspondingly reduced crj| s , as expected. However, since the gluon 
can only be calculated in the MS scheme for consistency we give our results 
in the MS scheme for both channels. 

The extremes we found that still allowed the cross section to be resummed 
are shown in the dotted curves of Fig. 3. Note that these calculations are 
not consistent with the results obtained in the previous studies and are only 
meant to indicate the possible range of cr res . The upper curves for cc produc- 
tion are obtained with m c = 1.3 GeV/c 2 , fi = m c and the GRV HO densities 
while the lower curves are calculated with m c = 1.8 GeV/c 2 , fi = 2m c and 
the MRS D— ' densities. For bb production, the upper curves are calculated 
with rrib = 4.5 GeV/c 2 , fi = vn^jl and the GRV LO densities. The lower 
curves are obtained with nib = 5 GeV/c 2 , \i = Im^ and the GRV HO densi- 
ties. (Since \x > Qq for the GRV HO and MRS D— ' the same scales are used 



6 We note that if the MRS D— ' distributions were refit with a smaller Qo so that /i = m c 
could be used, we would then expect the charm results to be similar to those with the 
GRV HO distributions. 



9 



in each, producing very similar results.) In Table 1 we show the appropri- 
ate value of a s for each parton density and scale used to calculate both our 
central results and our extreme cases. To illustrate how the resummation 
technique works at the extremes, in Figs. 4 and 5 we compare <r res and <r app 
for the upper and lower bounds given by the dotted curves in Fig. 3. 

In Fig. 4, the upper limit of cc production obtained using the GRV HO 
distributions produces convergence for /i = 0.25m c in the qq channel and 
Hq = 0.35m c in the gg channel. These values, with m c = 1.3 GeV/c 2 are not 
so different from those we found with m c = 1.5 GeV/c 2 . The smaller quark 
mass produces a stronger increase of <r res with energy, as suggested in Fig. 
3. The np — > cc cross section already diverges at \/S = 15 GeV in Fig. 3(a). 
The upper limit of bb production, calculated with the GRV LO distributions, 
converges for /i = 0.1 lm;, in the qq channel and 0.15m;, in the gg channel. 
Note that the smaller scale, used for GRV LO results in a smaller /i 

needed for resummation. However the variation of both cr rcs and cr app with 
fio is quite strong since at fi « 0.6m?,, so/m 2 « 1 and the separation of the 
S4 integral is no longer applicable. In this case, the bb cross section begins to 
diverge at large y/S, due also in part to the large value of a s for the lower 
scale. 

In Fig. 5, the lower limit of cc production obtained using the MRS D— ' 
distributions produces a convergent result for /i = 0.5m c in the qq channel 
and /io — 0.7m c in the gg channel. The larger value of m c makes the re- 
summation stable at higher energies. However, the resulting K^ p is so large 
this possibility appears to be ruled out. The lower limit of bb production, 
calculated with the GRV HO distributions, converges for // = 0.425m b in 
the qq channel and 0.625m;, in the gg channel. The ratio fio/m is similar for 
the lower limits of cc production since the scale is fi = 2m in both cases even 
though different parton densities were used. In general, it is the value of the 
scale that determines the hq at which the resummed and approximate results 
converge rather than the quark mass. For example, in bb production in pp 
interactions, the MRS D— ' and GRV HO distributions give similar results 
for fi Q since // = m& can be used for both p3| . Note that the higher mass and 
scale increases the relative qq contribution, particularly for ir~p interactions. 

To summarize, we have presented a comparison of the NLO and resummed 
cross sections for charm and bottom production using similar /io values to 
those obtained for top quark production. The resummation of the S+V 
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logarithms produces an enhancement of the NLO results. The resumma- 
tion technique appears to work reasonably well for m c = 1.5 GeV/c 2 and 
m b = 4.75 GeV/c 2 , resulting also in K^ p 1. At the moment, for a consis- 
tent description of both cc and bb production by pions and protons, the GRV 
HO distributions are the only choice available. We reach this conclusion con- 
cerning the GRV HO densities for two reasons: It is currently the only set 
with the same low x treatment of the pion and proton parton distribution 
functions. It also allows us to treat the scale on the same level for all heavy 
quarks, fi = ttiq, producing a uniform convergence of the resummed cross 
sections at the same Hq for charm, bottom, and top. Presumably another 
set of parton distributions that satisfies these criteria would produce results 
similar to those given in Fig. 2 and the central curves of Fig. 3. More precise 
data on the charm, bottom and top production cross sections in the thresh- 
old regions should clarify the situation and yield interesting information on 
the interplay between perturbative and nonperturbative physics. 
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cc 




GRV HO (2 loop) 


GRV HO (2 loop) 


MRS D-' (2 loop) 


A 3 (MeV/c) 


248 


248 


280 


m c (GeV/c 2 ) 


1.3 


1.5 


1.8 


Hj 'm c 


1 


1 


2 


a s 


0.3 


0.278 


0.204 


bb 




GRV LO (1 loop) 


GRV HO (2 loop) 


GRV HO (2 loop) 


A 4 (MeV/c) 


200 


200 


200 


m b (GeV/c 2 ) 


4.5 


4.75 


5.0 


fi/m b 


0.5 


1 


2 


a s 


0.312 


0.187 


0.155 



Table 1: The values of Af and a s for each of the parton densities and scales 
we consider. The left column gives the parameters for the upper bound on 
<7 res . The middle column shows the parameters used for our principle results. 
The right column displays the parameters that give the lower bound on cr rcs . 
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Figure Captions 

Fig. 1. Fractional contributions of the NLO channels to the total 0(a^) 
cc and bb production cross sections as functions of v^S"- We use m c = 1.5 
GeV/c 2 , m b = 4.75 GeV/c 2 , and the GRV HO parton densities with the two- 
loop corrected a s from PDFLIB fL2| . All the contributions are given in the 
MS scheme. We show the gg (solid), qq (dashed), and (q + q)g (dot-dashed), 
contributions to cc production in (a) n~p interactions and (b) pp interactions 
and bb production in (c) n~p and (d) pp interactions. 

Fig. 2. The fio dependence of the resummed cross section for cc and bb 
production. We use m c = 1.5 GeV/c 2 , m?, = 4.75 GeV/c 2 , and the GRV 
HO parton densities with the uncorrected two-loop a s . All the contributions 
are given in the MS scheme. We show cr T qq ? (solid), <r T g e g s (dashed), a^ p (dot- 
dashed) and a;* 13 (dotted) at = 15 GeV for cc production in (a) n~p and 
(b) pp interactions and at \/S = 30 GeV for bb production in (c) n~p and 
(d) pp interactions. 

Fig. 3. The resummed, improved, and NLO cross sections are shown as 
functions of v^S* for the cc and bb production. We show the total resummed 
cross section (solid) and the total improved cross section (dashed), both with 
fio = 0.15m for the qq channel and fio = 0.35m for the gg channel, and the 
total O(af) cross section (dot-dashed) results with /i = m using m c = 1.5 
GeV/c 2 and mj, = 4.75 GeV/c 2 . We also show the extreme values of <r res 
obtained when varying the quark mass, scale, and parton densities in the 
dotted lines. The results are given for cc production in (a) rr^p interactions 
and (b) pp interactions and bb production in (c) ir~p and (d) pp interactions. 

Fig. 4. The /^o dependence of the upper bound of the resummed cc and 
bb production cross sections shown in the upper dotted curves of Fig. 3. 
We show a v q f (solid), a v g c g s (dashed), a^ p (dot-dashed) and a^/ (dotted) at 

= 15 GeV for cc production in (a) ir~p and (b) pp interactions and at 
= 30 GeV for bb production in (c) ir~p and (d) pp interactions. For cc 
production, the curves correspond to the upper dotted curves in Fig. 3(a) 
and (b) with /i = m c = 1.3 GeV/c 2 and the GRV HO densities. For bb pro- 
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duction, the curves correspond to the upper dotted curves in Fig. 3(c) and 
(d) with 2 ix = m b = 4.5 GeV/c 2 and the GRV LO densities. 

Fig. 5. The /x dependence of the lower bound of the resummed cc and bb 
production cross sections shown in the lower dotted curves of Fig. 3. We show 
a T q f (solid), a r ™ (dashed), a^ p (dot-dashed) and (dotted) at y/S = 15 

GeV for cc production in (a) 7r~p and (b) pp interactions and at = 30 
GeV for bb production in (c) n~p and (d) pp interactions. For cc production, 
the curves correspond to the lower dotted curves in Fig. 3(a) and (b) with 
fi/2 = m c = 1.8 GeV/c 2 and the MRS D— ' densities. For bb production, 
the curves correspond to the lower dotted curves in Fig. 3(c) and (d) with 
fi/2 = m b = 5 GeV/c 2 and the GRV HO densities. 
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